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Abstract 

A free-energy functional for a crystal proposed by Singh and Singh (Europhys. Lett. 88, 16005 
(2009)) and which contains both the symmetry conserved and symmetry broken parts of the direct 
pair correlation function has been used to investigate the crystallization of a two-dimensional fluid. 
The results found for fluids interacting via the inverse power potential n(r) = e (a/r)" for n= 3, 6 
and 12 are in good agreement with experimental and simulation results. The contribution made 
by the symmetry broken part to the grand thermodynamic potential at the freezing point is found 
to increase with the softness of the potential. Our results explain why the Ramakrishnan-Yussouff 
(Phys. Rev. B 19, 2775 (1979)) free-energy functional gave good account of freezing transitions of 
hard-core potentials but failed for potentials that have soft core and/or attractive tail. 

PACS numbers: 64.70.D-, 63.20.dk, 05.70.Fh 
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I. INTRODUCTION 



Freezing is a basic phenomenon, the most inevitable of all phase changes. When a hquid 
freezes into a crystalline solid its continuous symmetry of translation and rotation is broken 
into one of the Bravais lattices. A crystalline solid has a discrete set of vectors Ri such 
that any function of position such as one particle density p(r) satisfies p(r + Ri) = p(r) 
for all Ri [1]. Because of localization of particles on lattice sites, a crystal is a system of 
extreme inhomogeneities where values of p(r) show extreme differences between its values on 
the lattice sites and in the interstitial regions. The density functional formahsm of classical 
statistical mechanics has been used to develop theories for liquid- solid transitions [2, 3]. 
This kind of approach was initiated in 1979 by Ramakrishnan and Yussouff [4]. 

The central quantity in this formulation is the excess reduced Helmholtz free energy due 
to interparticle interactions of both the crystal >lex[p] and the hquid Aex[pi\ [4,5]. For the 
crystal Aex[p] is a unique functional of p(r) whereas for the liquid Aex{pi) is simply a function 
of liquid density pi. The density functional formalism is used to write an expression for Ag^^fp] 
in terms of one and two-particle distribution functions of the solid [2-5]. The direct pair 
correlation function (DPCF) that appears in this expression is a functional of p(r) [2]. When 
this functional dependence is ignored and the DPCF is replaced by that of the coexisting 
uniform hquid [4,5] or by that of an " effective uniform fluid " [6,7] the free energy functional 
becomes approximate and fails to provide an accurate description of freezing transitions for 
a large class of intermolecular potentials. Attempts to include a term involving the three - 
body direct correlation function of the coexisting liquid in the free energy functional failed 
to improve the situation [8, 9]. 

It has recently been emphasized [10-12] that at the freezing point a qualitatively new 
contribution to the correlations in distribution of particles emerges due to spontaneous 
symmetry breaking. This fact has been used to write the DPCF of a frozen phase as a sum 
of two qualitatively different contributions; one that preserves the continuous symmetry of 
uniform liquid and the other that breaks it and vanishes in the liquid. The double functional 
integration in density space of a relation that connects Ae2:[p] to the DPCF led to an exact 
expression for Aej;[p]. The freezing transitions in three-dimensions have been investigated 
using this new free energy functional. The results found for the isotropic-nematic transition 
[10], crystaUization of power- law fluids [11] and the freezing of fluids of hard spheres into 
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crystalline and glassy phases [12] are very encouraging. 

In this paper we apply the free energy functional to investigate the liquid-solid transition 
in two-dimensions. It may, however, be noted that in contrast to three-dimensional solid, 
a two-dimensional solid melts in two-steps; the intermediate phase known as hexatic has 
a very narrow stability region in between liquid and crystal [13-15]. Since inclusion of the 
hexatic phase in the density functional formalism has not so for been possible, we neglect 
its presence and focus on the freezing of a fluid into the crystalline phase. Similar approach 
has been taken by others [16-21]. Here our motivation is to examine how well this new free 
energy functional (described briefly in the following section) compares with other free energy 
functionals in describing the crystallization of two-dimensional fluids. 

The paper is organized as follow : In Sec. II we give a brief description of the free-energy 
functional for a symmetry broken phase that contains both the symmetry-conserving and 
symmetry- broken parts of the DPCF. In Sec. Ill we describe methods to calculate these 
correlation functions for a two-dimensional system. The theory is applied in Sec. IV to 
investigate the freezing of power-law fluids into a crystalline solid of hexagonal lattice. The 
paper ends with a brief summary and perspectives in Sec. V. 

II. FREE ENERGY FUNCTIONAL 

The reduced free energy functional A\p\ of an inhomogeneous system is functional of p(r) 
and is written as [2] 



where A is a cube of the thermal wavelength associated with a molecule. The second 
functional derivative of the excess part Aex[p] with respect to p(r) deflnes the DPCF c(ri, 
of the system [2] , 



A\p] = Aid[p] + Ae:c\p]- 



(2.1) 



The ideal gas part Ai^ is exactly known and is written in terms of p(r) as 




(2.2) 



S^A,,[p] 



c(ri,r2, [p]). 



(2.3) 
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The function c that appears in tliis equation is related to the total correlation function 
/i(ri,r2) through the Ornstein -Zernike(OZ) equation, 

h{ri, rs) = c(ri, ra) + j dr3c{ri, r3)p(r3)/i(r2, rg). (2.4) 

Both functions h and c are functional of p(r) 

Since breaking of continuous symmetry of a uniform liquid at the freezing point gives rise 
to a qualitatively new contribution to correlations in the distribution of particles [10-12], 
the DPCF of the frozen phase is written as a sum of two contributions; 

c(ri,r2; [p]) = c^'Wr^ - v,\, po) + c^'\r,,r,; [p]), (2.5) 

where c*^°^ is symmetry-conserving and c^''^ symmetry-broken parts of the DPCF. While 
c^°^ depends on magnitude of interparticle separation r and is function of average density 
Po — {p{'^)), c^^^ is invariant only under discrete set of translations and rotations and is 
functional of p(r). 

Using Eq. (2.5) we rewrite Eq. (2.3) as 



Sp{ri)6p{r2) 



-c(°)(|r2-ri|,po), (2.6) 



5p{Yi)5p{r2) 

where A(^\,M + A^''\Ap\ = ^eApl The expressions for A^^^e^H and A^^'^exM are found from 
functional integrations of Eqs. (2.6) and (2.7), respectively. In this integration, as described 
elsewhere [10-12], the system is taken from some initial density to the final density along 
a path in the density space; the result is independent of the path of integration. These 
integrations give, 

A(°)e.[p] = AM + /3p - In (M) -\jdv^j dr2(p(ri) - pi){p{v2) - pM''\\v2 - ri|), (2.8) 
and 

A^'\Ap\ = -\ j dr, J dv2{p{v,) - po)(p(r2) - Po)c^'\r,, r^), (2.9) 
where 

c(°)(|r2 - ril) = 2 / 'dXX j 'rfA'c(o)(|r2 - ri|,p; + AA'(po - Pi)), (2.10) 
Jo Jo 



c^'\r,,r2)^4['dCc['dC' [ ' dXX [ ' dX' c^'\r,,r2; XX' polCC' Pa)- (2.11) 
Jo Jo Jo Jo 

Here Aex{pi) is excess reduced free energy of the coexisting uniform liquid of density pi 
and chemical potential fi , po = pi{l + Ap*) is average density of the solid and /3 = ksT, 
ks being the Boltzmann constant and T is the temperature and pc is an order parameter 
arising due to breaking of symmetry. 

The expression for Aea;*^°^[p] given by Eq. (2.8) is found from functional integrations 
when density pi of the coexisting fluid is taken as a reference. The expression for Aea;^^^[p] 
given by Eq. (2.9) is found by performing double functional integrations in the density 
space corresponding to the symmetry broken phase. The path of integration in this space is 
characterised by two parameters A and ^. These parameters vary from to 1. The parameter 
A raises the density from zero to final value p^ as it varies from to 1, whereas the parameter 
^ raises the order parameters from zero to their final values pa- The result is independent 
of the order of integration. 

The free energy functional for the symmetry broken phase is sum of Aj^^fp], Aea;^°^[p] and 
Aj''\p] Thus 

A\p\ = j dvp{v)[ln{p{r)K) - 1] + + /3(p - ln{piK)) j dv{p{v) - pi) 

-\ j dr, J drMn) - pi){p{r2) - pi)c^'\\r, - r,\) 

-IJdr, I dr2(p(r0 -po)(p(r2) -po)c(^n|ri,r2|) (2.12) 

where c^°^ and c^^^ are given, respectively by (2.10) and (2.11). In deriving Eq. (2.12) no 
approximation has been introduced. In the Ramakrishnan-Yussouff free energy functional 
c^^^ is neglected and c*^°^ is replaced by c^^\ 

In locating transition the grand thermodynamic potential defied as 

-W ^ A- I3p J drp{r), (2.13) 

is generally used as it ensures that the pressure and the chemical potential of the two 
phases remain equal at the transition. The transition point is determined by the condition 
Al^ — Wi — W — 0, where Wi is the grand thermodynamic potential of the coexisting 
liquid. Prom above expressions one gets the following expression for AW; 

AW^J dr[p{v) In ^ - (p(r) - pi)] -\ j dv,dv,{p{v,) - pi){p{v,) - pi)c^''\\v^ - v,\) 
-i|rfrirfr2(p(rO -po)(p(r2) -po)c(''n|ri,r2|). (2.14) 

5 



Minimization of AW with respect to p(r) subjected to the perfect crystal constraint leads 
to 



ln^ = 0+ / dn{p{r2)-pi)d^'\\T2-r,\)+ f dr^ipir,) - po)d^'\r,,r2), (2.15) 
Pi J J 



where 



and 



c(°)(|r2 - ril) = f'dXc^'Wr^ - v,\,pi + A(po - Pi)), 
Jo 

c^'\r,,r2) = Cdil^d\c^\v^,T2Apo,iPG). 
Jo Jo 

The value of Lagrange multiplier (f) in Eq. (2.15) is found from the condition 

One needs, in principle, the values of c*^°) and c^'^^ to calculate self consistently the value 
of p(r) that minimizes W. In practice however, one finds it convenient to do minimization 
with an assumed form of p(r). The ideal part is calculated using a form for p(r) which is a 
superposition of normalized Gaussians centred around the lattice sites; 



p(r) = -^exp[-a(r-R„)2], (2.17) 



where a is the localization parameter. For the interaction part it is convenient to use the 
Fourier expression, 

p(r) = po + E Pce^""" (2-18) 

where G are reciprocal lattice vectors (RLV's) of the lattice and pc — PoP'G are order 
parameters. Taking Fourier transform of Eq. (2.17) one finds p,G — e^~'^^/'^'^\ 

III. APPLICATION TO CRYSTALLIZATION OF POWER-LAW FLUIDS 
A. Potential model 

We consider model fluids interacting via inverse power pair potentials u{r) = e(o"/r)" ; 
where e, a and n are potential parameters and r is molecular separation. The parameter n 
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measures softness of the repulsion; n = oo corresponds to the hard disk and n = 1 to the 
one component plasma. Such repulsive potentials can be realised in colloidal suspensions. 
One such systems in two-dimensions has been provided by paramagnetic colloidal particles 
in a pendant water droplet, which are confined to the air- water interface [13]. By applying 
an external magnetic field perpendicular to the interface, a magnetic moment is introduced 
in the particles resulting into a tunable mutual dipolar repulsion between them. The pair 
interaction thus created is repulsive and proportional to r~^. The crystallization of this 
system has been investigated by van TeefTelen et.al. [20,21] using several versions of density 
functional theory (DFT). Other example where short range repulsion between molecules is 
found is microgel spheres whose diameter could be temperature tuned [15]. Most computer 
simulation studies on these systems suffer from the finite-size effects. In case of hard disks 
recently a large scale Monte- Carlo simulation, large enough to access the thermodynamic 
regime, has been performed [22]. The result confirms two-steps transitions from liquid to 
solid with the intermediate haxatic phase[23,24]. However, the liquid- hexatic transition, in 
contrast to the prediction of Kosterlitz-Thouless-Halperin-Nelson-Young(KTHNY) theory 
[23,24] is found to be first-order while the hexatic-sohd transition is second-order. The 
density functional theory predicts the liquid-solid transition to be first-order. 

In addition to being a pair potential that can be realized in a real system, it has a well 
known scaling property according to which the reduced excess thermodynamic properties 
depend on a single variable (or coupling constant) which for a two-dimensional system is 
defined as 

7 = (pa2)(/3e)2/" = p*r*-'/". 
Using this scaling the potential is written as 

where F = 7"/^ and r is measured in units of ao = {1/ pf^"^ . 

B. Calculation of c^'^^v) and its derivatives with respect to p 

The pair correlation functions of a classical system can be found in any spatial dimensions 
as a simultaneous solution of the OZ equation (Eq.(2.4)) and a closure relation that relates 
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functions h, c and the potential u(r). Several closure relations including the Percus- Yevick 
(PY) relation, the hypernetted chain (HNC) relation, modified hypernetted chain (MHNC) 
relation etc. have been used to describe structure of a uniform fluid [25]. We may however 
note that while the OZ equation is general and connects the total and direct pair correlation 
functions of liquids as well as of symmetry broken phases, the closure relations that exist in 
the literature have been derived assuming translational invariance [25]. They are therefore 
valid only for normal fluids. We use the integral equation theory involving suitable closure 
relations to find symmetry-conserving part of pair correlation functions h^^\r) and c^^\r) 
and their derivatives with respect to density. The symmetry-broken part of the DPCF is 
calculated using a method described in ref.[ll]. 

The OZ equation for a uniform system of density p reduces to 

/i(o)(r) c(°)(^) + pj c^r'c(°)(^>^°n|r' - r|). (3.1) 

The HNC closure relation and a closure relation proposed by Roger and Young [26] by 
mixing the PY and the HNC relations in such a way that at r = it reduces to the PY and 
for r ^ oo it reduces to the HNC relation, can be written together as 

exp{x(r)/(r)} - 1 



/i(o)(r) = exp{-/3u{r)) 



1 + 



1, (3.2) 



m 

where x(r) = /i^°)(r) — c*^°)(r) and /(r) — 1 — e~'^^ is a mixing function with an adjustable 
parameter < ^jJ < oo. For = oo or, f(r)= 1, Eq. (3.2) reduces to the HNC closure 
relation. In the Roger- Young relation, ip is chosen to guarantee thermodynamic consistency 
between the virial and compressibility routes to the equation of state. 

The differentiation of Eqs. (3.1) and (3.2) with respect to density p yields following 
two relations 

= Qp ^ j ^ )h^^\\Y -r\) + pj dr —^^h^^>{\r - r|) 
+p J dr c(°) (r) ^- ^ (3.3) 



and 



dp dp 

The closed set of coupled equations (3.1) — (3.4) have been solved for four unknowns h^^\ 
^(0) ^ ^^^^^(^) and The method can be extended to include higher order derivatives. 



In Fig. 1 we plot the Fourier transform of c^^^ (r) defied as 

c^°\q) = p J drc^^\r)e'''-\ (3.5) 

for (n, 7)= (3, 4.30), (6, 1.30) and (12, 0.90). The values given in Fig. 3.1 (a) for n=3 are 
in good agreement with values found by van Teeffelen et.al [20,21] (see Fig.l of their paper). 
As has been reported in ref. [21] the HNC closure underestimates values of c^^\q) whereas 
the Roger- Young (RY) closure gives relatively better but not very accurate values. In Fig. 
1 we also give values found from an approach proposed by Kang and Ree (KR) [27]. 
The exact closure relation which one finds from the liquid state theory [25] can be written 
as 

1 + /^(o) (r) = ^ W (r) = ex.p[-pu(r) +x{r) + B{r)], (3.6) 

where B{r) is the bridge function. In the HNC closure relation B{r) is taken equal to zero. 
In the KR approach the bridge function calculated for a reference potential and denoted as 
Bo{r) is used for B{r) in Eq. (3.6). The evaluation of -Bo(r) is done prior to and separated 
from the main integral equation by solving the Martynov-Sarkisov [28] integral equation. We 
briefly summarise here the way this is done for soft repulsive potentials in two-dimensions. 
The potential u{r) is first divided into a reference uo{r) and a perturbation part Up{r). 

u{r) = Uo{r) +Up{r), (3.7) 

where 

uo{r) — u{r) — F{r) if r < a 

= if r>a (3.8) 



Up{r) — F{r) if r < a 

= u{r) if r> a (3.9) 

Here F(r) — u{a) — u'{a){a — r) and a is nearest neighbour distance for hexagonal lattice at 
given density p. The Bo{r) for the reference potential is evaluated using the OZ equation 

ho^'\r) - co(°HO = Xo^'\r) ^pj dv co^''\r')ho^'\\v - r'|) (3.10) 
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and closure relation 



1 + /io^°^(r) = exp[-^iio(r) + xo(r) + So(r)]. 



(3.11) 



For -Bo(r) the Mortynov- Sarkisov [28] relation 



So(r) = [l + sxo(r)]'/^-l 



Xoir), 



(3.12) 



with s=2 is used. The values of Bo{r) are found by solving Eqs.(3.10)-(3.12) self-consistently. 
The value of -Bo(r) as a function of r for n= 3 is plotted in Fig. 2. The nature of -Bo(r) 
is same as was found in case of three-dimensions [27]. This value of -Bo(r) found for the 
reference potentials is used in relation (3.6) which is used to solve the OZ equation self- 
consistently to get values of h'^^^r) and c'^°)(r). The values of c^°^(g) found by this method 
are shown in Fig. 1 by full lines. These values are close to the simulation values given by 
van Teeffelen et.al [20,21] for n=3. For n=6 and 12 values found from the RY closure and 
values found from the KR method are close showing that for short-range repulsive potentials 
the RY closure yields good values of pair correlation functions. 

C. Calculation of c^'') (ri, ra) 

For a crystal c*^^)(ri,r2) is invariant only under discrete set of translations corresponding 
to lattice vectors {Rn}- If one chooses a center of mass variable Vc — ^'"^^'"'^ difference 
variable r = r2 — ri, the c^'') can be written as [11,12] 



where G are RLV's. Since c*^^^ is real and symmetric with respect to interchange of ri and 
r2 , c^'^\r) — c^~'^\r) and c^'^\r) — c^'^\—r). The function c*^^^(ri,r2) can be expanded in 
terms of higher body direct correlation functions of uniform liquid [2] ; 



+ -J rfr3rfr4C4^ ^(ri,r2,r3,r4;po)(p(r3) -Po)(p(r4) -po) H , (3.14) 

where p(r„) — pQ — Y.g Pae^^'^"^ , and c^^^^ are the n-body direct correlation functions 
of a uniform liquid of density po- These correlation functions are related to derivatives of 



c('')(rl,r2) = Eexp(^G.r,)c(^)(r;[p]), 



(3.13) 



G 
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■^^\r, po) with respect to density po as follows [2] ; 



92c(o)(r) 
dpo"^ 



j dr^ j rfr4C4(°)(ri,r2,r3, r4), (3.15) 



etc. 



The values of derivatives of 0*^°^ (r) appearing on the left hand side of above equations can 
be found using the integral equation theory described above. The usefulness of this method 
to find c*^^''(ri, r2) depends on convergence of series Eq. (3.14) which is a series in ascending 
powers of order parameters, and our ability of finding values of n-body {n > 3) direct 
correlation functions from Eq. (3.15). Barrat et.al [8] have shown that 03^°^ can be factored 
as C3(°)(ri, r2, r3) = t{ri2)'t{'i'iz)t{'i'2z) and the function t(r) can be determined from relation 



This method can be extended for higher Cn^°) [29]. Since c^^\ri, r2) is averaged over density 
p and over order parameters pa the contributions made by successive terms of Eq. (3.14) 
in 74ea;*^*''[p] is expected to decrease rapidly [11]. In the case of three dimensions it was 
found that it is only first term of the series (3.14) which needs to be considered to describe 
accurately the fluid -solid transition [11, 29]. In two-dimensions the convergence is expected 
to be faster as number of nearest neighbours is less compared to the three-dimensions and 
therefore the higher body correlation functions expected to be less important. In view of 
this, we consider here the first term of the series (3.14) only and examine its contribution in 
stabilizing the hexagonal lattice at the transition point. 

From known values of ^^-p — we solved numerically Eq. (3.16) to find values of t(r) for 
different values of 7. In Fig. 3, we plot values of t(r) for n = 3, 6 and 12 at values of 7 close 
to the freezing point. 

Taking only the first term of Eq. (3.14) the expression for c^^\v-i, Y2) in terms of t(r) can be 



(seeEq. (3.15)) 




(3.16) 



written as 



c(^)(ri,r2) =poE / dv^t{\v^-v,\)e'''-^H{\v^-V2\). 



(3.17) 



G 
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Using the relation 

1 1 

ra = ^{ri + ra) + (rg - ri) - ^{^2 - ri) 

= rc + r --(r2-r), 
Eq. (3.17) reduces to Eq. (3.13), i.e. 

cW(ri,r2) = ^e'^"'=c(^)(r), 

G 

where 

c(G)(r) = pol^Gt{r)e-'^'^-' j dv t{r')e'^-'' t{\Y - r|). (3.18) 

Using the relation e^^ '" = Em(«)'"'^m(Gr)e'"'('^G-0r) ^^ere Jm{Gr) is the Bessel function of 
the first kind of integral order m we find following expression for c^^\r)] 

c(^)(r) = ^(^)*^CM^''^^)e"''^^^ (3-19) 

M 

where 

CM^''\r) ^ PoPGt{r)Y,Brn{r,G)Jrn+M{\Gr)e''''^^ (3.20) 

and 

S„(r, G) = y dkkt{k)Jm{kr)e''^'^^ J dr'r' Jm{kr')Jm{Gr') (3.21) 

For hexagonal lattice M = 0, ±6. The value of CM^^\r) depends on values of order param- 
eters pg and on the values of RLV's. In Figs. 4-6, we plot harmonic coefficients c^'^'^^r) 
and CQ-^^r) for n = 3, 6 and 12 for RLV's of first four sets, respectively. For different set 
of RLV's cm^'^\t) varies with r in different way; the values in all cases become negligible 
for r ( measured in units of ao = (l/p)^^^)> 1.5 . For any given value of G, the values of 
CQ^^\r) is about an order of magnitude larger than CQ^'^\r) at their maxima and minima. 
As magnitude of G increases value of CM''^\r) decreases and after the sixth set of RLV's 
values of CM'^'^\r) become negligible for all values of n. 

IV. LIQUID- SOLID TRANSITION 

Substituting expressions of p(r) given by Eqs. (2.17) and (2.18) and of c(^)(ri, ra) given 
by Eq. (3.19) in Eq. (2.14) we find. 

N N ^ N ^ N ' ^ ' 



12 



where 




(4.2) 



N 

AWq 
N 



1 ' ^1 

-o 51 51 /^Gi/^-G-GiC(Gi + -G). 



-lA7c('^)(o)-iEl/^Gre^°HG), 



(4.3) 



(4.4) 



N 



G Gi 



where A7 = ( 2:^-^ j ; the subscripts s and 1 represent sohd and hquid respectively. 



Hear AWj^, A Wo and AW^ are respectively, the ideal, symmetry-conserving and 
symmetry-broken contributions to AW. The prime on a summation in Eq (4.4) indicates 
the condition G 7^ 0, d 7^ and G d 7^ 0, and 



We used above expressions to locate the liquid-crystal (hexagonal lattice) transition by 
varying values of 7, A7 and a. The results given in Table 1 for n= 3, 6 and 12 correspond 
to the RY closure relation. We note that the contribution arising due to symmetry broken 
part of the DPCF is far from neghgible and its importance increases with the softness of 
potential. While it is about 7.3 % to the symmetry conserving term for n=12, it increases 
about 44 % for n=3. This explains why the Ramakrishnan-Yussouff theory gives good 
results for hard core potentials but fails for potentials that have soft core and/or attractive 
tail. As the contribution of is negative, it stabilizes the sohd phase. Without it the 
theory strongly overestimates the stabihty of fluid phase specially for softer potentials [20,21] . 
The contribution made by the symmetry broken part of the DPCF is, as expected, small 
compered to that in three-dimensions (3D) at the freezing point for the same potential. For 
example, the contribution in 3D [11] for n=12 is 22.2 % compared to 7.3 % in 2D whereas 
for n=6 the contribution is 37 % in 3D and 18 % in 2D . 

In Table 2 we compare results of the present calculation using both the RY closure and 
the KR procedure to calculate pair correlation functions for n=3 with the results found from 
other free-energy functionals as reported in ref [21]. The experimental results obtained from 
real-space microscopy measurements of magnetic colloids confined to an air-water interface 
[13] and values found from numerical simulations [30,31] are also given in the table. While 
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the RY closure gives slightly higher values of Fj and Tg compared to experimental values, 
the KR closure gives slightly lower values. But these values along with the values of other 
parameters, particularly the values of AF — Fg — Ff, are in better agreement with the 
experimental values compared to any other versions of the DPT. Although the extended 
modified weighted-density approximation (EMA) [32] with Verlet closure [33] gives values 
oiFf which is close to the one found by us using the KR closure, but the values of AF is 
significantly lower; AF ema = 0.16 compared to the values found by us AF = 0.41, and the 
experimental value, 0.75 . 

The real-space experimental data are not available for other systems. The computer 
simulation results [34, 35] show the liquid-solid transition at 7; = 1.51 and 0.986 respectively 
for n = 6 and 12. These values are close to the one given in Table 1. 

V. SUMMERY AND PERSPECTIVES 

We used a free energy functional that contains both the symmetry-conserving part of 
the DPCF c(o)(r) and the symmetry-broken part c(^)(ri,r2) to investigate the freezing of 
a two-dimensional fluid into a two-dimensional crystal of hexagonal lattice. The values of 
c^°^(r) and its derivatives with respect to density p as a function of interparticle separation 
r have been determined using an integral equation theory comprising the OZ equation and 
the closure relations of Roger and Young [26] and of Kang and Ree [27] . For soft potential 
(n=3) the two results are found to differ; the KR closure seems to give better result. For 
more repulsive potentials the two results are close as shown in Fig. 1. For c*^''^(ri, which 
is functional of p(r) and is invariant only under discrete set of translations and rotations, 
we used an expansion in ascending powers of order parameters. This expansion involves 
higher body direct correlation functions of isotropic phase, which in turn were found from 
the density derivatives of c^°)(r) using a method proposed by Barrat et.al [8]. 

The contribution of symmetry-broken part of DPCF to the free energy is found to depend 
on nature of pair potentials; the contribution increases with softness of potentials. This result 
is in agreement with that found in three-dimensions and explains why the Ramakrishnan- 
Yussouff free-energy functional was found to give a reasonably good description of the freez- 
ing transition of hard core potentials but failed for potentials that have soft core and/or 
attractive tail. The results found here and the results reported for 3D indicate that the 
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theory described here can be used to investigate the freezing transitions of all kinds of 
fluids. 

Since our free energy functional takes into account the spontaneous symmetry breaking 
it can be used to study various phenomena of ordered phases. The results indicate that 
the density-functional approach provides an effective frame work for theoretical study of a 
large variety of problems involving inhomogeneities. However, the question not adequately 
addressed yet is the size of fluctuations effect which play important role in two-dimensional 
systems. The other important question is the inclusion of hexatic phase in the theory. 
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Tabic 1. Freezing parameters q;,7;,A7 and the pressure P at coexistence along with 
the contributions of ideal, symmetry-conserving and symmetry-broken parts of These 
results correspond to the Roger- Young closure [26] . 



n 


a 


li 


A7 


N 


AWo 

N 


AWt 

N 


p 


3 


100 


4.96 


0.025 


2.50 


-1.74 


- 0.76 


75 


C 


100 


1.55 


0.010 


2.50 


-2.10 


- 0. 10 


31 


12 


96 


1.00 


0.050 


2.49 


-2.32 


- 0.17 


22 



Table 2. Freezing parameters r/(= 7/"''^), ^s{— 7s "''^) and the width of coexistence region 
AF — Tg—Vf, and the relative displacement parameter 2/q;) at the coexistence obtained 
from various density functional schemes. The MWDA stands for modified weighted density 
approximation, EMA for extended modified weighted density approximation, RY and KR 
refer to, respectively, the Roger- Young closure [26] and the Kang and Ree [27] closure. 
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FIG. 1: The dimensionless Fourier transform c^^\q) of the direct pair correlation function of c^^\r) 
plotted against kao (clq = {l/p)^^"^^ shown are data found from the integral equation theory using 
the RY closure (dashed line), HNC closure (dotted line) and KR closure (full line) at values of 7 
shown in figure (a), (b) and (c) for n=3, 6 and 12, respectively. 
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FIG. 3: Function t(r) vs r for n=3, 6 and 12 at the values of 7 above the fluid-solid transition 
point, r is in units of oq = (l/p)^^"^. 
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FIG. 4: Harmonic coefficients CM^^\r) for RLV's of first four sets for n=3, 7 = 4.42(r = 9.32). 
Notations are as follow: Full line represents values of the first set, the dotted line of the second set 
, dashed line of the third set and dashed-dotted line of the fourth set. The distance r is expressed 
in unit oq, where ao(= 
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FIG. 5: Harmonic coefficients CM^^\r) for RLV's of first four sets for n=6, 7 = 1.52(r = 3.52). 
Notations are same as in Fig 4. 
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FIG. 6: Harmonic coefficients CM^^\r) for RLV's of first four sets for n=12, 7 = 1.05(r = 1.34). 

Notations are same as in Fig 4. 
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